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ABSTRACT 

Aims. We further characterize the structures tentatively identified on thermal and chemical grounds as the sites of dissipation of tur- 
bulence in molecular clouds (Papers I and II). 

Methods. Our study is based on two-point statistics of line centroid velocities (CV), computed from three large '^CO maps of two 
fields. We build the probability density functions (PDF) of the CO line centroid velocity increments (CVI) over lags varying by an 
order of magnitude. Structure functions of the line CV are computed up to the 6'* order. We compare these statistical properties in two 
translucent parsec-scale fields embedded in different large-scale environments, one far from virial balance and the other virialized. 
We also address their scale dependence in the former, more turbulent, field. 

Results. The statistical properties of the line CV bear the three signatures of intermittency in a turbulent velocity field: (1) the non- 
Gaussian tails in the CVI PDF grow as the lag decreases, (2) the departure from Kolmogorov scaling of the high-order structure 
functions is more pronounced in the more turbulent field, (3) the positions contributing to the CVI PDF tails delineate narrow fila- 
mentary structures (thickness ~ 0.02 pc), uncorrelated to dense gas structures and spatially coherent with thicker ones (~ 0.18 pc) 
observed on larger scales. We show that the largest CVI trace sharp variations of the extreme CO linewings and that they actually cap- 
ture properties of the underlying velocity field, uncontaminated by density fluctuations. The confrontation with theoretical predictions 
leads us to identify these small-scale filamentary structures with extrema of velocity-shears. We estimate that viscous dissipation at 
the 0.02 pc-scale in these structures is up to 10 times higher than average, consistent with their being associated with gas warmer than 
the bulk. Last, their average direction is parallel (or close) to that of the local magnetic field projection. 

Conclusions. Turbulence in these translucent fields exhibits the statistical and structural signatures of small-scale and inertial-range 
intermittency. The more turbulent field on the 30 pc-scale is also the more intermittent on small scales. The small-scale intermittent 
structures coincide with those formerly identified as sites of enhanced dissipation. They are organized into parsec-scale coherent 
structures, coupling a broad range of scales. 

Key words. ISM: clouds, ISM: magnetic fields, ISM: kinematics and dynamics, turbulence 

1. Introduction in laboratory experiments of incompressible turbulence, in the 

atmosphere, and in the solar wind (see re cent reviews by 
Star formation proceeds via gravitational instability in dense lAnselmet et alJl200lt iBruno & Carboi^l2005l) . It manifests it- 
gas, but the respective roles of turbulence and magnetic ^elf i) as non-Gaussian tails in the probabiUty distribution of 
fields in that process are still debated issues fCiolek & Basu| ^jj quantities involving velocity differences (e.g. gradient, shear, 
2006t iTassis & Mouschovias 2004; Mac Low & Klesseii||2004t vorticity, rates of strain, and energy dissipation rate), ii) anoma- 



| Padoanetal. || 2001|; |Klessen| | 2001(: |Bate et al.||2002,) in spite of j^^^ g^^jj^g the high-order s tructure functions of the veloc- 



dedicated observational studies of magnetic fields in molecular increments (lAnselmet et al J flfjl . a nd Hi) the existence of 

clouds ( [Matthews & Wilson 2000; Crutcher 1999) and extensive coher ent structures of intense vorticity (I Vincent & Meneguzzil 

theoretical and numerical works devoted to characterizing the Hf^lMoisv & Jimenez '2004. hereafter MJ04). It has long been 

properties of the turbulence ( |Boldyrev et a l. 2002; Padoan et al. | unclear whether these three manifestations, which refer either to 

12003). A hybrid paradigm is taking shape, where turbulence j^e statistical properties of the flow or to its coherent structures, 

dominates the diffuse ISM dynam ics and niagnet ic fields gain vvere related 
importance as the scale decreases (ICrutcheiil20()5l) . Turbulence 

and magnetic fields are recognized as powerful stabilizing agents Statistical models make theoretical predictions in terms of 

in molecular clouds, and a critical issue remains, turbulence dis- the two-point statistics of the velocity field, with no link to any 

sipation: where, when, and at which rate and scale does it occur ? physical structures in the turbulent flow. In particular, the struc- 

A generic property of turbulent flows is the space-time ture functions S ,,(1) - {[v(r + I) - v(r)]'') of the velocity field 

intermittency of the velocity field. Intermittency is observed are expected to be power laws S ,,(1) oc /f*'''. As p increases, the 

structure functions give more weight to intense and rare events. 

Send offprint requests to: e-mail: philybla(aobs . uj f-grenoble . fr In principle then, a detailed description of the velocity field could 

* Based on observations carried out with the IRAM-30m telescope, be achieved with the knowledge of all Sp for p ^ oo, however, 

IRAM is supported by INSU-CNRS/MPG/IGN. the number of points Np needed to compute S p grows with p. 
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Fig. 1. Reprocessed IRAS map of the Polaris Flare 
jMiville-Deschenes & Lagachd l2005h . The map size is about 
10° X 10° or 27 X 27 pc. The parsec-scale field analyzed in this 
paper is shown as a white rectangle. The 100//m, 60/im and 
12yL/m emissions are red, green and blue respectively. 



and an educated guess is A^^ ~ 10''^^. Fortunately, theoretical 
models show that, with a limited number of orders (p > 3), some 
properties of the turbulence can still be tested. The Kolmogorov 
theory of no n-intermittent turbulence (hereafter K41, see e.g. 
iFrischI [r995') predicts ((p) - p/3, while experiments show a 
clear departure from this scaling, generally with ^(p) < p/3 for 
p > 3. The departure from the p/3 scaling is usually interpreted 
as a definition of intermittency. Statistical models predict inter- 
mittent scalings ({p) + p/3. One such model (She & Leveque 
11994 hereafter SL94, see Appendix A3) has an intermittent 
scaUng = p/9 + 2[1 - (2/3 )''^^1, in exce llent agreement with 
tunnel-flow experimental data (iBenzi et al.lll993h . 

The structural approach is inspired from laborat ory ex- 
perim ents showing filaments of high vorticity (Doua dv et al.l 
[l99l . Localized regions of high vorticity and rate of strain 
(and thus energy dissipation) are found in numerical simula- 
tions at high resolutio n of both incompressible and compress- 
ible (Porter etal. 1994) and magneto-hydrodynamical (MHD ) 
turbulence ( Vincent & Meneguzzi 1991; Mininni et al. 2006bl) . 
Recently, MJ04 found that intense structures of vorticity and 
rate of strain are respectively filaments and ribbons that are not 
randomly distributed in space but that instead form clusters of 
inertial-range extent, implying a large-scale organization of the 
small-scale intermittent structures. In 1024^' numerical simula- 
tions of inc ompressible turbulence , with variable large-scale stir- 
ring forces. iMininni et al.l (l2006al) have shown that the large and 
small-scale properties of the flow are correlated, namely that i) 
more intense small-scale gradients and vortex tubes are concen- 
trated in the regions where the large-scale shears are the largest, 
and that ii) the departure from the Kolmogorov scaling is more 
pronounced in these regions. They infer from these results that 
the statistical signatures of intermittency are linked to the exis- 
tence of the small-scale vortex tubes. 

Investigations of interstellar turbulence are plagued by pro- 
jection eff'ects that affect our knowledge of the velocity. Direct 
inversion of the observations being a highly degenerate proce- 
dure, progress relies on astrophysical observables derived from 
numerical simulations of tur bulence, and their confron tation to 
real data (see the review of lElmegreen & Scaloll2004l) . Along 
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Fig. 2. Map of the cold dust emission in the Taurus molec- 
ular complex built fr om the rep rocessed 60/vm and 100/vm 
IRAS maps (i Hilv-Blan t et al.l l2007b) following the method of 
lAbergel et al.l (119951) . The orientation of the projection of the 
magnetic fields on the sky (lieiles 2000 ) is shown as white 
segments of length proportional to the polarization degree. The 
parsec-scale field discussed in this paper is centered at aim) - 
04''40'"08.84^ (52000 = 24°12'48.40" and is shown as a white 
rectangle at offsets (0°,-l°). 



these lines, iLis et al.l d 19961) have shown that it is possible to 
trace the projection of the most intense velocity-shears with a 
measurable quantity based on two-point statistics of the velocity 
field: the line centroid velocity increments (CVI). This method, 
applied to different fie lds observed in CO lin es, a star-forming 
regio n dLis et al.ll 19981) and quiescent regions dPetv & Falgarond 
I2003L hereafter PF03) revealed filamentary structures unrelated 
with the gas mass distribution. The line centroid velocities (and 
their increments) are sensitive to optical depth effects and to 
density, temperature, and abundance fluctuations along the line 
of sight. Their relevance in any analysis of the statistical prop- 
erties of the actual velocity fields have therefore been repeat- 
edly questioned. Several recent studies have clarified the is- 
sue (Lazarian & Pogosyan 200 0; Miville-Deschenes et al. 20031 
ILevr ier 2004: lOssenkopf et al J 12006) without having provided 
any final answer to that question yet. The present work is part of 
a broad study in which we characterize the structures of largest 
CVI on thermal, chemical, and statistical grounds, and then re- 
peat this analysis in different turbulent clouds. The goal of this 
broad program is to shed light on what these structures actually 
trace. 

One of the two fields studied in this paper is the parsec- 
scale environment of two low-mass de nse cores in the Polaris 
Flare. In the vicinity of the dense cores, iFalgarone et al.l (l2006h 
(hereafter Paper I) find HCO^ abundances locally far in ex- 
cess of steady-state chemical predictions. The large measured 
abundances are consistent with a scenario that involves an im- 
pulsive heating of the gas, during which a warm chemistry is 
triggered, followed by a slow chemical and thermal relaxation. 
According to this scenario, the observed abundances and gas 
densities, in the range 200 and 600 cm"^, correspond to a gas 
that has cooled down to 100-200 K. In Hily-Blant & Fal garong 
( l2007ah (hereafter Paper II), we analyze the mass distribution 
of the gas in the environment of the dense cores and disclose 
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Fig. 3. Histograms of the equivalent width (Aveq - 

j r(v.,) dvjrpeak) of the i2cO(l - 0) line in Polaris (shaded 
area) and Taurus. The result of the log-normal fit to the Polaris 
histogram is also shown. 



localized regions where the '^CO profiles exhibit broad wings. 
Multi-line analysis shows that gas in these linewings is optically 
thin in the '^CO(l - 0) line and warmer than 25 K with density 
< lOOOcm"^. The new result is that this warm gas component 
is not widespread but localized in filaments. The dispersion of 
the orientation of the filaments of dense gas in the field, com- 
pared to the velocity dispersion, suggests that the turbulence 
in t hat field is trans- Alfve nic. In the forthcoming and last pa- 
per dFalgarone et al.ll2008l) . we will report '^00(1 - 0) IRAM- 
Plateau de Bure observations performed in this field, revealing 
milliparsec-scale structures with large velocity-shears. 

The present paper is dedicated to the statistical and struc- 
tural analysis of turbulence in similar parsec-scale translucent 
samples of gas, which belong to two different large-scale en- 
vironments: one, far from virial balance and devoid of stars, is 
the above-mentioned field in the Polaris Flare; the other lies at 
the edge of the Taurus-Auriga molecular complex. The statistical 
analysis is based on two-point statistics of the CO line emission 
observed at high spectral resolution, and the structural analysis 
consists in characterizing the regions of largest Une-centroid ve- 
locity increments. 

After presenting the data in Section 2, the characteristics of 
the turbulence in these two fields are derived based on the PDF 
of line CVI and structure functions (Section 3). We then show, in 
Section 4, that the regions of largest CVI on small-scales are fila- 
ments uncoiTelated with the distribution of matter, but coiTelated 
with the filaments of gas optically thin in '^CO(l -0) where large 
HCO^ abundances are found. It is shown that these filaments re- 
main coherent at the parsec-scale. In Section 5, we show that the 
ensemble of results, based either on the statistical or structural 
approach, provides a consistent description of the intermittency 
of turbulence in these two fields. In Section 6, the comparison 
of the radiative cooling of these structures with the fraction of 
the turbulent energy susceptible to being dissipated there further 
supports the proposition that the filaments of largest CVI some- 
how trace extrema of velocity-shear in the fields and pinpoint the 
sites of intermittent dissipation of turbulence. 



Table 1. Dispersions crgc (in kms ') of the PDF of CVI com- 
puted from the '^CO(l - 0) transition in the Polaris and Taurus 
fields for different lags I (in pixels). 
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Polaris 


Tauras 


[pixels] 


[pc] 


kms"' 


kms"' 


3 


0.02 


0.11 


0.05 


6 


0.05 


0.17 


0.07 


9 


0.07 


0.22 


0.09 


12 


0.09 


0.25 


0.10 


15 


0.11 


0.27 


0.11 


18 


0.14 


0.29 


0.12 




0.14 


0.30 






0.27 


0.46 




54(«) 


0.41 


0.57 




72(") 


0.54 


0.64 




90(«) 


0.68 


0.69 




108*°' 


0.81 


0.72 





(a) Computed from the KOSMA ^^CO{2 - 1) data of lBensch et"!!] I20Q1I) . where one 
pixel conesponds to 6 pixels of the IRAM data. The adopted distance is 150 pc for 
both fields. 



2. Observations 

2.1. Polaris field 

Paper II describes the observations and data reduction of the 
IRAM '-CO and i^^CO(l - 0) and (2 - 1) maps of the molecular 
cloud MCLD 123.5+24.9 located in the Polaris Flare. The loca- 
tion of the field mapped is shown in Fig. [T] The fully-sampled 
w 1 pc X pc maps cover the non-star-forming translucent envi- 
ronrnent of two dense c ores dGerin et al.lll997t iFalgarone et alJ 
119981; lHeithausenll2002l) . corresponding to » 3000 independent 
spectra (i.e. spaced by one beamsize or 20" at 115 GHz). The 
spectral resolution is 0.055km s"'. Since the signal-to-noise ra- 
tio of the J - 2 - I data is insufficient, only the J - I - 
transition are used in the present paper. For comparison with 
larger scale properties, we use the fully-sampled ^^CO(2 - 1) 
data from Bensch et al. (2001) obtained at a lower angular reso- 
lution (HPBW = 120") with the 3m KOSMA antenna. 

2.2. Taurus field 

The second field is located at the edge of the Taurus-Auriga 
molecular complex (Fig. |2]i. Observations were done with the 
IRAM-30m telescope. Observational strategy and data reduction 
are similar to those of the Polaris field. The maps, centered at 
(ffMOO = 04*'40™08.84^ 52000 = 24°12'48.40"), are fully sam- 
pled in '^CO and '^CO(l - 0). The data will be presente d in 
more detail in a later paper (iHilv-Blant & Falgaronell2008l) but 
here we give the properties relevant to the present work. A total 
of 1200 independent spectra was obtained with the same spectral 
resolution of 0.055 km s"' . In a small region (around 0,0" offsets 
in the map of Fig.|9]l, spectra show two separate components (at 
vlsr ~ 5 and 10 km s"'), and the analysis presented in this paper 
focuses on the main component at vlsr ~ 5 kms"', by blanking 
out the area corresponding to the high-velocity component. 

2.3. Comparison of the two fields 

Both fields are translucent. The visual extinction lies between 
Av - 0.6 and 0.8 mag at a resolution of 8' in the Polaris field 
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dCambresv et al.ll2001h and between Ay = 1 and 1.2 ma g in 
the Taurus field at the same resolution dCambresvl 1199 9*). As 
shown in Fig. [3] the most probable equivalent width (Aveq - 
J T(vt) dVf /Tpeak, X being the coordinate along the line of sight) 
of the '^CO(l - 0) line is a factor two larger in Polaris than in 
Taurus. However, since the lines are stronger in the Taurus field, 
the integrated intensities in both fields are similar. Not only is the 
equivalent width larger in the Polaris field than in the Taurus one, 
but so is the dispersion of these equivalent widths. This factor 2 
between the equivalent width translates into a factor 4 in the spe- 
cific kinetic energy ratio between the two fields. The distribution 
of the equivalent width in Polaris is also very well-fitted by a log- 
normal distribution centered at Avgq = 2 kms ' with dispersion 
1.2 kms~^ The parsec-scale velocity gradients, deduced from 
the centroid velocity maps, in the Polaris field 2 km s ' pc ') 
is also twice larger than in the Taurus field. 

These two parsec-scale fields are similar with respect to their 
size and average column density but have specific kinetic en- 
ergies that differ by a factor 4. Moreover, they belong to two 
very different environments on the scale of ~ 30 pc, that of 
Figs. [1] and |2] The total gas mass M,o, - 4.4 x lO'* Mq at 
the scale of 30 pc is close to the virial m ass in the Taurus- 
Auriga field dUngerechts & ThaddeusI [T987h . while it is more 
than six times lower in the Polaris Flare, M,„, - 5500 Mq with 
My = 3.6 X 10"* Mo dHeithausen & ThaddeusI [T990h . It is in- 
teresting that the virial masses of the two large-scale fields are 
close, because their velocity dispersion on the scale of 30 pc are 
similar, 3.8 and 4.8 km s"', respectively. 

In summary, the less turbulent parsec-scale field lies on the 
far outskirts of the virialized Taurus-Auriga molecular complex, 
while the more turbulent field belongs to a much less massive 
complex, far from virial balance. 



3. Two-point statistics of the centroid velocity 

3. 1 . Probability density functions of the line centroid velocity 
increments 



Following iLis et aP dl996l) . we analyze the two-point statistics 
of the centroid of the line-of-sight projection of the velocity V j, 
which we note as C{y, z) - C(r), where (y, z) is the position on 
the sky: 



C(r) = j T(r,v,)v,dvJ j r(r,v,)dv. 



(1) 



Increments of the centroid velocity between 2 points separated 
by / are defined by SC{r, I) = C(r + I) - C(r). This quantity will 
be called centroid velocity increment (CVI). The main difficulty 
of this method concerns the computation of C, which not only 
depends on the bounds of the integrals in Eq. [1] but is also af- 
fected by the signal-to-noise ratio (SN R). To circumvent the bias 
introduced by spatial noise variations. iRosolowskv et al] d 19991) 
degrade all the spectra to a unique threshold SNR. A different 
approach (PF03) has been adopted here, which uses the SNR of 
the integrated area as the optimization criterion to determine the 
spectral window used to compute C. The reason is that we have 
checked that the noise is already homogeneous in the data cubes, 
mostly as a result of the observing strategy consisting in several 
coverages of individual sub-maps in perpendicular directions. 

For a given value of I -\l\, we compute maps of CVI for each 
direction ///. A probability density function (PDF) is built from 
these maps by normalizing the histogram of CVI to a unit area. 
We thus obtain a PDF for each / and each direction. For each /, 





-5 5 -5 5 -5 5 

c5C, / C7(|5CJ 

Fig. 4. Normalized PDF, P„{SCi), of the centroid velocity in- 
crements (CVI) computed from the '^CO(l - 0) IRAM map in 
Polaris . The PDF are computed for different lags between pairs 
of points: Z = 3 to 25 pixels, and normalized to unity dispersion, 
such that the x-axis is in units of the rms o^c for each distri- 
bution. Only the bins containing more than 10 data points were 
kept. The values of the dispersion crgc are given in Table [1] A 
Gaussian of unit dispersion is also shown (dotted curve). 
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Fig. 5. Same as Fig.|4]for the '^CO(l -0) data towards the Taurus 
field. 



a PDF is computed with the CVI from all directions I /I, which 
we denote as P{6Ci). In order to get PDF with zero average and 
unit standard deviation and to ease the comparisons, we use the 
normalized PDF PniSCi). All bins of the 'P„{6Ci) associated to 
a number of points less than a given value Nj„i„ are blanked (see 
Appendix A.l). In the following, all P„(6Ci) have 32 bins, and 
the adopted minimum number of data points for a bin to be sig- 
nificant is A^min - 10. In a second step, for a given I, we compute 
the azimuthal average of the absolute value of the CVI, resulting 
in a single CVI map. In practice the structures seen in the non- 
averaged maps are not smeared out, though they appear thinner 
in some cases. 

Figs. |4] and |5] show the P„{SCi) computed for various lags 
from Z = 3 to 25 pixels in the Polaris and Taurus fields, respec- 
tively. The lag / = 3 is the shortest distance between two inde- 
pendent points (since the sampling is half the beam size), and 
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Fig. 6. Same as Fig . |4]for the 
lBenschetal.l(l200li 



^CO(2 - 1) KOSMA data from 



Table 2. Standard deviations cr (in km s of the centroid veloc- 
ity PDF ipos) and of the average Hne profiles (los), computed in 
three ways in the two fields (crj, 0-2, 0-3). 



Field 


Type 


Size 


0-1 


0-2 


0-3 


{^) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Polaris ' 


pos 


2.1 X 2.8 


0.54 


0.57 


0.60 


0.57 + 0.03 




los 




1.13 


1.10 


1.16 


1.13 ±0.03 


Polaris^ 


pos 


0.7 X 0.6 


0.31 


0.32 


0.25 


0.29 ± 0.04 




los 




0.97 


1.02 


1.10 


1.03 ±0.05 


Taurus'^ 


pos 


0.4 X 0.7 


0.12 


0.13 


0.11 


0.12 ±0.01 




los 




0.50 


0.60 


0.66 


0.59 ± 0.07 



(1) 

(2) 
(3) 
(4) 
(5) 
(6) 
(7) 



For the Polaris field, the IRAM (small) and KOSMA (large) datasets are taken sepa- 
rately, t: based on the '-C0(2 - 1) data, ^^C0{\ - 0) data 
the type of PDF {pos or los) 

Map sizes (in pc X pc) are computed assuming a distance of 150 pc 

a"i is the standard deviation 

0-2 is derived from a Gaussian fit 

o"3 = A\'|.q/2.35 

average of the three determinations 



I - 25 corresponds to the largest lag with significant number of 
pairs of points. The number of data points corresponding to the 
three most extreme bins for I - 3 and 25 are in the range 12-50 
and 20 - 500, respectively. 

The PDF at large lags (/ > 15) (Fig.|4|i are nearly Gaussian 
and become slightly asymmetrical at / - 12, an effect we at- 
tribute to large-scale velocity gradients. Such effects cancel out 
at lags smaller than the characteristic scale of these gradients, 
hence the more symmetrical shape of the PDF at small lags. It 
is not obvious that such large-scale gradients should be removed 
(see discussion in PF03). The fields mapped here are expected 
to be small with respect to the integral scale of turbulence L, at 
least on the order of the molecular cloud size itself. As the lag de- 
creases from / = 25 to Z = 3 pixels, non-Gaussian tails develop. 
These tails are more pronounced in the Polaris field than in the 
Taurus one, with CVI values up to 6 times the dispersions ctsc of 
the unnormalized PDF (see Table[T]i. However, since the number 
of points in Taurus is lower than in the Polaris field, the mini- 
mum level of probability reached is an order of magnitude higher 
(10^^ instead of 10"^ in Polaris). We also computed the PDF of 
the increments for the large-scale KOSMA data in the Polaris 
field (Fig. |6]l and we also find increasing non-Gaussian tails as 
the lag decreases. The dispersions of the PDF are reported in 
Table [1] and are seen to smoothly connect with the small-scale 
values computed in the IRAM field. 

In both fields, the dependence of cr^c with / can be well- 
fitted by a power law ctsc °^ f''^- The dispersions crgc in the 
Taurus field are a factor ^ 2 smaller than in the Polaris field. 
This ratio is also found when comparing the velocity dispersions 
- either across the plane of the sky {pos) or along the line of sight 
{los) - in the two fields (see Table |2|i. Furthermore, the ratio of 
the los to pos dispersions suggests that the depth of the cloud 
is larger than the extension in the plane of the sky {los > pos) 
jOssenkopf & Mac Lowir2002l) . 

3.2. Non-Gaussianity: the flatness 

The deviations of the PDF from a Gaussian shape can be quan- 
tified using the flatness (or kurtosis) of a distribution, defined by 



where the pth-order moment (for even p) is computed as 
<5Cf) = J SC^ rjdCi) d{6C,). The flatness equals 3 for a 
Gaussian distribution. The uncertainties, here and for all the mo- 
ments computed in what follows, are estimated by using two 
'P„{6Ci) with two thresholds, A^min = 1 and 10, and by computing 
the mean and rms of the two outputs. Fig.|7]displays the flatness 
of the PDF at all lags, for the two fields. For the Polaris field, the 
IRAM and KOSMA data have a flatness close to 3 at large lags, 
which confirms the visual Gaussian shape of the coiTesponding 
PDF in Fig. |4] The flatness increases at smaller lags as a result 
of non-Gaussian tails. For the Taurus field, however, the flatness 
stays close to 3, confirming that the non-Gaussian tails are less 
pronounced. 

3.3. Structure functions of the line centroid velocities 

By analogy w ith studies performed on the velocity field {e.g. 
IShe et al.ll2001 ). we computed the structure functions of the line 
CV, using the PDF of the centroid velocity increments P„{6Ci), 
a procedure that allows a check of the convergence of the struc- 
ture functions by filtering out doubtful points in the PDF. The 
structure functions are evaluated by a direct integration of the 
PDF: 



5p(0= mfPMSCil) d{\6C,\). 



T{1)^ 



(2) 



(3) 



Structure functions of velocity are frequently normalized to 
the third-order function for two reasons: first because in incom- 
pressible, isotropic, and homogeneous turbulence, ^(3) = 1 is 
an exact result; second because oscillations in S p{l) - I plots are 
damped when the S p are plott ed against 5 3, a p roperty called ex- 
tended self-similarity (ESS) (Ben zi et al.ll993h . so that the range 
of scales over which the structure functions are power laws is 
wider. 

Calculations of high-order structure functions are suscepti- 
ble to errors since, as p increases, any spurious large fluctuation 
largely affects S p. The calculation of S p based on Eq. fallows 
us to reject bins populated by too small a number of points, and 
to study the influence of iiTelevant bins of the PDF (for which 
N < A^min = 10). We also applied the procedure described in 
i Leveque & Shd d 19971) and Padoan et al. (2003) for determining 
the highest significant order: the peak of the histogram of {dC'i'l 
occurs for a value of 6C1 that must be represented by a significant 
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Table 3. Exponents ^j, - ^(p)/^(3) of the ESS structure func- 
tions of the CV for the Polaris and Taurus fields (see Fig.[T4ll. 



1 



10^ 
i [pixels] 



Fig. 7. Flatness of the CVI (see Eq.|2]i computed from the nor- 
maHzed PDF of CVI. Squares: Polaris '^00(1 - 0) IRAM data. 
Triangles: 1^00(2- l)KOSMA data. Circles: Taurus 12CO(1-0) 
IRAM data. 



number of points in the PDF of 6Ci. The maximum significant 
order found is p = 6. 

The scalings ofSp with 53 are shown in Fig.[8]for the Polaris 
and Taurus fields. They are power laws. Exponents of the struc- 
ture functions are calculated by fitting the ESS diagrams for lags 
wider than 2 pixels and smaller than 30 to 60 (see Fig.[8]l. Error 
bars on the exponents (1 - 3%) are estimated from weighted 
averages of the results from two calculations coiTesponding to 
A^tnin=l and 10. The ESS exponent values are given in Table [3] 
The absolute values of the second-order structure function in the 
Taurus field are a factor 6 lower than in Polaris, confirming that 
the Taurus field contains less kinetic energy. The non-ESS fits 
of the CV structure functions lead to values for ^(3) = 1.35 and 
1.60 in Polaris and Taurus, which differ significantly from the 
expected value (^(3) = 1) for the velocity in incompressible, 
homogeneous, and isotropic turbulence tKjitsuk et al 2007.) . 



4. The spatial distribution of tlie largest line 
centroid velocity increments 

In the following section, we discuss the spatial structures of the 
largest CVI in the different maps. We compare maps of CVI 
computed on large and small-scales with large and small beams. 
For the sake of simplicity, we call shear the CVI value divided 
by the lag over which it is measured, 6C1/I. This will be justified 
at the end of Section 5. 



4.1. Locus of the extrema of CVI in the IRAM fields 

Figs. |9] and IA.3I show the maps of azimuthally averaged CVI 
({\6Ci\}) computed for two lags (1 = 3 and 18 pixels) in both 
fields. The grey scale is the magnitude of the azimuthally aver- 
aged CVI at a given position on the sky. Because of the averaging 
procedure, the exact values of the CVI in those maps are not triv- 
ially related to the values of the non-averaged PDF of Figs. |4]and 
12] However, regions of large CVI in the maps do correspond to 
the positions responsible for the non-Gaussian tails in the PDF 
of Figs. 141161 
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(a) From the large scale '^CO(2 - 1) data from Bensch etalj|2001l 



(b) SL94 and B02 are the lShe & Levequel jl994i) and iBoldvrev et alj 120021) scalings of 
the velocity structure functions (see Section 5.2). 



At a small lag (/ = 3), in both fields, the spatial distribution of 
the bright regions with large CVI delineates elongated structures. 
When the lag is larger (18 pixels), the contrast of the structures 
above the background values fades away. Yet, in the Polaris field, 
the structure around (-1000", -200") is still visible with / = 18, 
and for the two lags of 3 and 18 pixels, the largest CVI is located 
in the northwestern corner of the map. 

These maps show that the positions of the largest CVI, in 
Polaris and Taurus, are not randomly distributed but are con- 
nected and form elongated structures. In the Taurus field, the di- 
rection of the most prominent CVI structure is parallel to the 
projected orientatio n of magneti c fields measured in the NE 
corner of the field (lHeilesll2000l) (see also Fig. [B- In Polaris, 
the scatter of their orientations relative to the magnetic field is 
larger (Paper II). In both Taurus and Polaris, their character- 
istic half-maximum width, measured on transverse cuts, is re- 
solved (30" after deconvolution from the lag, or d ^ 0.02 pc), 
and their aspect ratio is often greater than 5. The surface frac- 
tions covered by the regions where the {\6Ci\) are larger than 
3craver (whcrc cTaver IS the dispcrsion of the distribution {\6Ci\}), 
are 10% and 28 % in Polaris and Taurus, respectively. Last, the 
non-averaged CVI in these structures (see Section 3. 1) are 5 and 
4 times larger than the dispersion osc (see Figs. |4| and |5]l in 
Polaris and Taurus, respectively. For / = 3, the coiTesponding 
shears are 5 x 0.11 kms"V0-02 pc = 30 kms"' pc"' in Polaris 
and 4 x 0.05 kms"V0.02 pc = 10 kms"' pc"' in Taurus. The 
most turbulent field on the parsec-scale (Polaris) is therefore that 
where the largest small-scale shears are measured. 

4.2. Comparison of the extrema of CVI with the CO emission 

We stress here that our work is based on the statistics of the 
extrema of CVI (called E-C VI in what follows), unlike what is 
done in most analyses (e.g. lEsquivel & Lazarianll2005l) . where 
the full distribution of CVI is considered. We show below which 
specific features of the CO line emission are associated with the 
extrema of CVI. 

Fig. [TOj displays two '^CO(l - 0) space-velocity cuts made 
across the Polaris map at longitude offsets -500 and -800". The 
run of the averaged CV and CVI along these cuts is shown to 
illustrate that the largest CVI are mostly due to very localized 
broad CO linewings. As expected, however, some of the large 
variations in the centroid velocities due to these broad wings are 
reduced by opposite variations due to fluctuations in the line- 
core emission. 

This is seen better in Fig. (TTj where the extrema of cen- 
troid velocity increments in the Polaris field are overplotted on 
the '^CO wing emission and the '^CO integrated emission. The 
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Fig. 8. Structure functions S p(l) plotted against ^3(0 (p - 1, . . . , 6) for Polaris (left for the IRAM data and middle for the KOSMA 
data) and Taurus (right). A power law is fitted to each order for I > 2 pixels, and for I < 60, / < 30, and Z < 40 from left to right. 



'^CO wing emission (top panel) is optically thin and associ- 
ated with warm and tenuous gas emission (Paper II). The '^CO- 
integrated emission (bottom panel) is used here as a proxy for the 
molecular gas column density. While the largest CVI are not spa- 
tially correlated with the '^CO-integrated emission, they closely 
follow the boundaries of the optically-thin '^CO emission in the 
broad '^CO linewings. 

Given the high latitude of the Polaris cloud, the structures 
responsible for the broad '^CO linewings most likely belong to 
the Polaris Flare (Fig. [TJ. The non-Gaussian tails of the PDF 
at small lags are thus associated to local structures on the scale 
of 30 pc. Their distance, and therefore their size, is known to 
within 20%, while the lags of the PDF decrease by an order of 
magnitude (from 25 to 3). The / - 3 pixel PDF is thus sensitive 
to true small-scale structures created by the turbulence in the 
Polaris Flare. This comparison shows that the regions of largest 
shear are not associated with the bulk of the condensed matter 
in the field, but are instead correlated with warmer and more di- 
luted gas. We therefore infer that, unlike centroid velocities in- 
crements in general, the E-CVI we analyze are not due to density 
fluctuations or radiative transfer effects in optically thick gas. 

We now address the issue of the chance coincidence of un- 
related pieces of gas on the line of sight. In the two translucent 
fields, CO is not expected to trace the full molecular content of 
the clouds, essentially as a result of photodissociation processes . 
This has possibly been observed bv ' Sakamoto & Sunadal (l2003h 
who show discontinuous CO emission in low extinction regions 
in the Taurus complex. These spots of CO emission are how- 
ever embedded in the underlying turbulent molecular gas, un- 
detected because mostly made of molecular hydrogen, and pre- 
sumably continuous. The velocity field deduced from the CO 
emission lines thus carries the statistical properties of that turbu- 
lent molecular gas. 

Nonetheless, projection effects are inevitable, and a key pa- 
rameter is the ratio //L of the lag / over which CVI are measured 
to the unknown depth of the cloud along the line-of-sight L. In 
their work on 5 1 2'' num erical simulations of mildly compress- 
ible turbulence. Lis et al.l(ll996l) have computed PDF of CVI for 
a lag of 3 pixels, for which this ratio is 3/512=0.006. They show 
that the E-CVI trace extrema of {(V xv)y}^^^ + {(V xv),}"^^^. Since 
this quantity is a los integration of a signed quantity (the two 
projections of the vorticity in the plane of the sky), its extrema 
are due to a fe w exceptional va lues present on the line of sight. 
For this reason iLis etal.l([T996h say that the E-CVI trace the pro- 



jection of large velocity-shears (or vorticity) in turbulence. In 
Polaris and Taurus, L is not known but we conservatively adopt 
a value in the range 1-30 pc. For the smallest lag (/ = 3 pixels) 
the ratio is in the range l/L = 0.001 - 0. 02. Our observa tional 
study thus falls into the regime tested by Lis et al.l (11996 ). and 
the filaments associated with the E-CVI thus trace the projection 
of regions of extreme velocity-shear somewhere on the line of 
sight. The observed value of the shear, though, is of course an 
upper limit. We are therefore confident that E-CVI trace genuine 
extrema of line-of-sight velocity fluctuations. 

4.3. Parsec-scale coherence of the regions of largest CVI: 
IRAM and KOSMA data 

We compare here the properties of the E-CVI (values and spa- 
tial distribution) from the KOSMA and IRAM data sets of the 
Polaris field. The comparison is not, however, straightforward 
for two reasons. First, the value of the centroid velocity is af- 
fected by the beam size, and second, the computation of the CVI 
filters out any structures that are much larger than the lag. 

Fig.[T2]displays the spatial distribution of the CVI computed 
on large scales with the KOSMA data, for a lag 1=3 pixels 
(180"). Regions of large increments are spatially resolved fila- 
ments: cuts across the structures provide an average thickness of 
200" deconvolved from the lag, or 0.18 pc. These structures are 
about 7 times thicker than those found in the IRAM field. The 
prominent KOSMA structure around (123.29°, 25.1 1°) smoothly 
connects with the northwestern IRAM structure (contours from 
Fig.fTTI). This is seen more clearly in Fig.[T3]where the values of 
the non-averaged CVI along this structure are displayed. Fig.[T3] 
illustrates three points: /) the CVI from the KOSMA and IRAM 
datasets decrease monotonously from north to south along this 
structure, //) the IRAM CVI measured over / = 180" are all larger 
than those measured with KOSMA over the same physical lag, 
and the CVI measured with the KOSMA telescope over a lag 
of 3 pixels (180") are similar to those measured at the same po- 
sitions with the IRAM telescope with a lag 6 times smaller (30") 
and . 

The latter property is unexpected: it suggests that the CVI 
are similar in this region whether they are measured with a small 
beam and a small lag (IRAM) or a large beam and a large physi- 
cal lag. If the KOSMA structures were only due to beam-dilution 
of the IRAM ones, the KOSMA CVI for I = 180" would be 



8 



Hily-Blant P. et al: Intermittency in diffuse molecular gas 




Fig.9. Spatial distribution of the azimuthally averaged CVI (i.e. (|5C/|)) computed with lags/ - 3 based on the '^CO(l-O) line. The 
grey-scale gives the (|i5C/|) in kms"'. The dark regions correspond to large values of the CVI associated to the tails of the P{6Ci) 
(Fig.|4]and|5]l. The orientations of the magnetic fields (Heiles 2000) are also shown. Left panel: Polaris field. Contours indicate the 
0. 11 and 0.22 km s"' levels. The largest CVI (0.30 km s"') appear in the NW corner of the field. The 4 crosses indicate the positions 
where the HCO^ abundances have been measured (Paper I). Right panel: Taurus field. Map center is a'2000 - 04''40™08.84% (Saooo = 
24°12'48.40". Offsets are in ". Contours: 0.05 kms"' with 0.01 kms"' steps. Note that the largest CVI (0.10 kms"') in that field 
are 3 times smaller than in the Polaris field. The blanked areas around (0, 0) offsets correspond to the positions where the 10 km s"^ 
component is present. 



smaller than the IRAM ones for / = 30". In other words, the 
KOSMA CVI structures are real and are sub-structured: the same 
velocity variations (< 0.5 kms"' for positive offsets) are mea- 
sured on small (IRAM) and large (KOSMA) scales. 

The largest velocity-shears at the KOSMA resolution are 
5x0.30/0.18 ~ 8.3kms"' pc"'. The surface fraction covered by 
the large increment structures where (|(5C;|) > 3craver are close, 
10 and 16% for the IRAM and KOSMA data, respectively. These 
two fields, observed with different telescopes and resolutions, 
thus show similar statistical properties and demonstrate the co- 
herence at the parsec-scale of the structures of largest centroid 
velocity increments. This is discussed in the next section. 

5. The two facets of intermittency: statistical and 
structural 

In Section 3, the two-point statistics of the line centroid ve- 
locities were found to display marked non-Gaussian behaviors. 
In Section 4, the emission responsible for these non-Gaussian 
statistics is resolved into coherent structures. We here compare 
these statistical and structural characteristics with theoretical 
predictions and recent numerical results regarding the intermit- 
tency of turbulence. 

5.1. Self-similarity of the centroid velocity increments in the 
Polaris field 

The IRAM and KOSMA PDF from Figs.|4]and|6]bear an appar- 
ent contradiction: the PDF built with the IRAM data with a lag of 
150" (15 pixels) is nearly Gaussian, while the KOSMA PDF with 
a similar physical lag of 1 80" (3 pixels) is not: non-Gaussian tails 



in the KOSMA PDF originate from the largest CVI structures 
like the most prominent one discussed in Section 4.3. In Fig. [131 
we show that the bulk of the IRAM CVI for / = 180" (dark 
dots) are below 3crsc - 0.9km s"'. It is the particular location 
of the IRAM field with respect to the CVI maxima seen in the 
KOSMA field that prevents the detection of a number of CVI in 
excess of 3crgc large enough to depart from Gaussian statistics. 
Would the IRAM field be centered closer to the CVI maximum 
in the KOSMA data (around 123.29°, 25.11°), a larger number 
of CVI in excess of 3o-sc might have been measured. 

However, in both fields, the non-Gaussian tails of the PDF 
of CVI grow as the lag decreases. This behavior is routinely 
observed in laboratory and numerical experiments, where it is 
interpreted as a signature of the intermittency of the velocity 
field. In such experiments, effects of the non-Gaussian statistics 
are visible with either the transverse or the longitudinal veloc- 
ity increments (Frisch 1995; Mininni et al. 2006a). The present 
analysis shows various degrees of non-Gaussianity: the statistics 
in Taurus are nearly Gaussian (flatness close to 3), while both 
Polaris data sets show clear departure from Gaussianity. If, fol- 
lowing PF03, we attribute the non-Gaussian statistics of the CVI 
to the intermittency of the turbulence in the two sampled molecu- 
lar clouds (see below), the new result, here, is that intermittency 
is as pronounced at a lag / - 180" in the large field (KOSMA 
PDF) as it is at lag / = 30" in the IRAM field. In both datasets, 
the non-Gaussian tails extend to 5.5 - 6crsc. This suggests that 
intermittency is not only a small-scale phenomenon but that it 
is present and has the same statistical properties on a scale six 
times larger. 

As mentioned in the introduction, similar conclusions have 
been reached by MJ04, who find that the intense structures of 
vorticity and rate of strain in hydrodynamical turbulence form 
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Fig. 10. Position-velocity cuts at constant {u offsets (A^u - 
-500", left and Ain = -800", right). Grey-scale: main-beam 
intensity in K. Dashed curve: centroid velocity. Full curve: az- 
imuthally averaged CVI for a lag I = 3 pixels (with an additional 
offset of -6.5 km s"' for clarity). 



clusters of inertial range extent, implying a large-scale organiza- 
tion of the small-scale intermittent structures. 



5.2. The intensity of small-scale intermittency versus 
large-scale shear 

We now compare the scaling of the /^'''-order structure functions 
of CV with p in the framework of the SL94 model. Structure 
functions are increasingly sensitive to the tails of the CVI PDF 
(E-CVI) as the order p increases. Since we have shown (Sect 4.2) 
that the E-CVI stem from velocity fluctuations, it is interesting 
to compare the scaling of high order structure functions of the 
CV to theoretical predictions based on the velocity field. 

The SL94 model has three parameters (see Appendix A. 3). 
One of the three parameters, < /? < 1, describes the level of 
intermittency: /3 —> I corresponds to the non-i ntermittent cas- 
cade with - pi 3. The two other parameters dBoldvrev et al.l 
I2OO2I) describe the scalings of velocity in the cascade and the 
dimension D of the most intermittent structures. In the SL94 
model, the scaling of the velocity is v/ ~ Z'^^. It assumes that 
the most intermittent structures are filaments (D - \) and that 
the level of intermittency is jS = 2/3. The associated ESS struc- 
ture function exponents 1p - ^{p)/({3) are then predicted to be 
^f- = p/9 -\-2[l - (2/3)''/^]. According to this class of models, 
as the level of intermittency increases, the ESS exponents be- 
come smaller than p/3 for p > 3 and the departure from the K41 
scaling increases with p. Following the SL94 approach, further 
theoretical models were developed for compressible and magne- 
tized turbulence (Politano & Pouquet 1995; Miiller & Biskamp 
and tested against numerical simulations. IBoldvrev et aXI 
I propose a similar scaling to describe compressible MHD 
turbulence assuming sheet-like intermittent structures (D = 2) 
and a more intermittent cascade (/? = 1/3) but do not allow 
dissipation of large-scale modes in shocks. They find ^^'^^ - 
pl9-\-\ - (1/3)''^"' (hereafter called the B02 scaling), in excellent 
agreement with numerical simulations of super-Alfvenic MHD 
turbulence. However, this scaling has never been tested against 
observations of the turbulent velocity field of molecular clouds. 

In Fig. [141 we compare the scalings of the CV structure func- 
tions exponents (see Table O with the predictions of the SL94 
and B02 models for the velocity field. In the Polaris field, for 
either data set, the measured exponents follow the SL94 scaling 
closely but differ significantly from that of B02. This apparent 
agreement with the SL94 scaling is unexpected since this model 
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Fig. 11. Contours of the CVI (1 km s ') computed in the Polaris 
field for 1-3 pixels (see Fig. |9ll overplotted on the '^CO(l - 
0) emission integrated in the velocity range [-2,5 : 0] kms"' 
i.e. optically thin '^CO emission (top) and the '^CO integrated 
intensity tracing the bulk of the dense gas (see Paper II) (bottom). 



describes incompressible and unmagnetized turbulence, every- 
thing interstellar turbulence is not. However, the effect of line- 
of-sight averaging in the CV structure functions is not known. 
Further interpretation of the underlying physics requires con- 
frontation with CV structure functions based on numerical sim- 
ulations. The exponents in the Taurus field do not follow any of 
the three scalings, and are halfway between the non-intermittent 
K41 scaling {((p) - p/3) and the SL94 scaling. The Taurus field 
is thus less intermittent than the Polaris one, a result consistent 
with the flatness measure of the non-Gaussianity of the 'P{6Ci). 
This result is also in fine with the recent numerical findings of 
iMininni et al. I (I2006al) . who show that the characteristics of the 
large-scale flow play an important role in the development of 
small-scale intermittency and determine its statistical properties. 

The use of CV structure functions to determine the three pa- 
rameters of the SL94 class of models not only requires a large 
number of data points but also "calibration" of the weighting 
performed by C V upon the velocity field, using numerical simu- 
lations of compressible turbulence. The values that we have de- 
termined from the exponents ({p) may be useful, though, and we 
give them in the Appendix. 

In summary, the line CV exhibit the statistical and structural 
properties characterizing the intermittency of the velocity field 
in theoretical models and numerical simulations of turbulence: 
1) the non-Gaussian statistics of the CVI at small lag, 2) the 
self-similarity of structures of largest CVI and the existence of 
inertial-range intermittency, 3) the anomalous scaling of their 
high-order structure functions similar to wthat is found for the 
structure functions of velocity. Last, we find that the more in- 
termittent field on small scales has the larger dispersion of non- 
thermal velocity on large scales. 

The above properties are borne by the non-Gaussian tail of 
the CVI PDF, which we have shown to be associated with pure 
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Fig. 12. Spatial distribution of CVI (km s~') c omputed on large 
scales using the data of .Bensch et al. I (l200l ( grey-scale with 
levels indicated on the right scale, dotted contours indicate the 
0.1 kms ' level). Full-line contours indicate the CVI computed 
on small scales in the lRAM-30m data cube (see Fig.fTTT). Note 
the spatial coincidence of the IRAM structures at the tip of the 
prominent KOSMA structure at (123.4°, 25°), where the CVI are 
0.30kms-i. 



velocity fluctuations. They support our proposition that statis- 
tics of the '^CO line centroid velocity and, more specifically, 
their E-CVI may be used to disclose the statistical and structural 
properties of intermittency in the underlying velocity field. In 
what follows, we therefore ascribe the E-CVl to the intermittent 
structures of intense shears. 



6. Discussion 

6. 1 . Influence of gravity 

To test the role of gravit y in the g enerat ion of non-Gaussian 
statistics of the velocity, iKlessenI (l2000h built the two-point 
statistics of the velocity field, in SPH numerical simulations of 
turbulence, both with and without self-gravity . The author com- 
pared the numerical PDF of CVI with observed PDF and con- 
cluded that the inclusion of self-gravity leads to better agree- 
ment with the observed PDF in molecular clouds. It was fur- 
ther argued in PF03 that this result was indeed expected since 
the observed regions used in the comparison are forming stars, 
hence the importance of self-gravity. The situation is drastically 
changed in the two translucent fields we have studied: they do 
not form stars, are far from any such regions, and both are non- 
self-gravitating on the parsec-scale of our observations. Indeed, 
we have shown that the field with the more prominent non- 
Gaussian tails, namely Polaris, is located at high latitude and 
is embedded in a larger structure (the Flare) far from virial bal- 
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Fig. 13. Values of the CVI in the Polaris field, along the NW- 
SE CVI structure from (123.2, 25.2 °) to (123.6, 24.9 °) seen in 
Figin] CVI ai-e from KOSMA data for I = 3 (diamonds), IRAM 
/ - 3 (open circles), and IRAM / = 18 (filled circles). Offsets 
increase from NW to SE, with the zero position corresponding 
to the NW corner of the IRAM field (see Fig. lA.2l ). 



ance. This strongly supports that, in the type of fields we analyze, 
gravity is not at the origin of the PDF tails. 

Nonetheless, gravity is the ultimate source of gas motions in 
the universe, from galaxy clusters to GMC and collapsing cores, 
so it cannot be ignored. If the cloud mass were distributed in 
tiny cloudlets of very high density that would rarely collide, then 
gravity might play a significant role in the gas velocity statis- 
tics. Here, we assume that the fluid approximation is valid, and 
because the Reynolds number is so large, the gas motions are 
turbulent, by definition. 

6.2. The intermittency of turbulence dissipation 

In numerica l experiments such as those of MJ04 (see also 
ISreenivasan| [T999), maxima of energy dissipation are found, at 
small-scale, in the vicinity of the vorticity filaments. A large 
fraction of the dissipation of turbulence may be concentrated 
in the regions of largest CVI, the small-scale intermittent struc- 
tures. We illustrate this point with estimates of energy transfer 
based on our observations in the Polaris field. 

The transfer rate in the cascade on scale / is e; = l/2pvf/Z 
with V/ the characteristic turbulent velocity fluctuations on that 
scale. This assumes that the time to transfer energy from scale I 
to smaller scales writes as r/ - l/vi. At the parsec-scale of the 
Polaris cloud 

fL = 5.5xlO"^Vipc erg s"' cm"^ (4) 

where the average density p = jjNn/L and the velocity disper- 
sion vl = 1.5 kms ' are those derived in Paper II. There, Lpc is 
the unknown depth along the line of sight, expressed in pc. The 
CO cooling rate averaged over the whole field on the parsec- 
scale is (Paper II) 

Aco = 10-24/Lpe ergs-'cm-l (5) 

These two rates are thus close, because the depth along the 
line of sight is not much greater than 1 pc (see Section 3.1), 
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Fig. 14. Values of the exponents of the structure functions (see Fig.[8]l with error bars (see Table[3]). For comparison, the K41 (full), 
SL94 (dashed), and B02 (dotted) scaUngs are indicated. From left to right; Polaris (IRAM), Polaris (KOSMA), Taurus. 



and may lead to the conclusion that CO is able to radiate the 
turbule nt energy away, confirming the results of IShore et aP 
( 1200 6'). Actually, this would be true if the cascade were fill- 
ing space uniformly or, in other words, for a non-intermittent 
cascade. From the previous section, this assumption is cer- 
tainly not valid. Indeed, the cumulative distribution of 6C^ for 
I - 3 (see Fig lB.lt shows that the points with CVI larger than 
3o-sc represent only 2.5% of the total, while they contribute to 
25% of the total of 6Cf. The velocity field has two contribu- 
tions (solenoidal and dilatati onal) to the energy d issipation rate 
e = -J- ( |Vxv p + ■ vp) (Kritsuk et al."2007'). As shown by 
iLis et ai.l d 19961) and iPetv & Falgarone (2000) based on numer- 
ical simulations, the E-CVI serve as a proxy for large vorticity 
regions. Thus, assuming that the energy dissipation on scale I is 
proportional to 6Cf, the cumulative distribution of Fig. lB.ll sug- 
gests that the local dissipation rate at Z = 3 (or ~ 0.02 pc) in the 
E-CVI regions, is akeady 10 times larger than the average rate 
over the field, fE-cvi > lOe/.. Note that these numbers are about 
the same for the large-scale Polaris field, while in the less inter- 
mittent Taurus field, these E-CVI represent only 1 % of the total, 
still contributing to 5% of the dissipation. 

On the actual dissipation scale, presumably smaller than 
0.02 pc, the local dissipation rate is still higher by an unknown 
factor. Since the turbulent dissipation is not space-filling, it in- 
duces high local heating rates. This suggests that other cool- 
ing agents, e.g. the pure rotational lines of H2 or the fine struc- 
ture line of C^, may be dominating the cooling in these regions 
(Paper I andl Falgarone e t al.ll2007h. Observations are still lack- 
ing that would allow a comparison of the turbulent transfer rate 
with the CO cooling rate on scales smaller than 0.02 pc. 

In spite of the self-similarity of the intermittent structures 
discussed in Section 5, the bulk of the dissipation is likely to 
take place in the smallest structures. The largest CVI are propor- 
tional to (Tgc and thus to The corresponding shears there- 
fore scale as r^^^ providing an observed scaling of the dissipa- 
tion rate with lengthscale I, e/ oc r\ Now, we use the finding 
of MJ04, who show that the tails of the probability distribution 
functions of the volume of individual dissipative structures (ei- 
ther intense vorticity or strain-rate) decrease approximately as 
p(V) ~ V^^. Whether these structures are cylinders {V oc f ) 
or sheets (V oc I), the integrated dissipation is therefore always 
dominated by the dissipation which takes place on the smallest 
scales, because p(V)ei oc in the first case or oc in the sec- 
ond case. 



This confirms the important point for the evolution of molec- 
ular clouds that dissipation of turbulence is concentrated in a 
small subset of space. The induced radiative cooling, and there- 
fore the dissipation rate, have to be searched on scales on the or- 
der of the milliparsec in emission lines more powerful than the 
low-y CO transitions. The value of the rate itself may thus be 
directly observable in line emissions (pure rotational lines of H2, 
C^) that can only be distinguished from UV-excited emission by 
observations at very high angular resolution. 

7. Conclusion and perspectives 

We performed a statistical analysis of the turbulence towards two 
translucent molecular clouds based on the two-point probability 
density functions of the '^CO(l - 0) line centroid velocity. 

Thanks to the excellent quality of the data, we prove the non- 
Gaussian tails in the PDF of the line centroid velocity increments 
on small-scales, down to a probability level of 10 We show 
that the largest CVI, in both fields, delineate elongated narrow 
structures (~ 0.02 pc) that are, in one case, parallel to the local 
direction of the magnetic field. In the Polaris field, these fila- 
ments are well-correlated to the warm gas traced by the optically 
thin '^CO(l -0), while they do not follow the distribution of mat- 
ter traced by the '■'CO. Using large-scale data, we have shown 
that these filamentary structures remain coherent over more than 
a parsec. Furthermore, the similar statistics found in the IRAM 
and KOSMA maps of this field suggest that both samples be- 
long to the self-similar turbulent cascade. In the Polaris field, the 
high-order structure function exponents, computed up to order 
p - 6, significantly depart from their Kolmogorov value. 

Through the properties of the tails of their PDF, i.e. the E- 
CVI, the line centroid velocities in these two clouds are found 
to carry the main signatures of intermittency borne by a turbu- 
lent velocity field. The departure from the Gaussian statistics of 
the centroid velocity increments on small-scales is therefore as- 
cribed to the intermittency of turbulence, i.e. the non-space fill- 
ing character of the turbulent cascade. The structures of largest 
CVI trace the intermittent structures of intense shears and the 
sites of intermittent turbulence dissipation. We show that these 
intermittent structures, on the 0.02 pc-scale, harbour 25% and 
5% of the total energy dissipation, in the Polaris and Taurus 
fields, respectively, although they fill less than 2.5% and 1% of 
the cloud area. We find that both fields are intermittent and that 
the more intermittent velocity field on small scales (the Polaris 
field) belongs to a molecular cloud far from virial balance on 
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the scale of 30 pc. In contrast, the less intermittent (the Taurus 
field) belongs to a virialized complex. The more turbulent field is 
thus the more intermittent. Interestingly enough, the less turbu- 
lent field is embedded in a star-forming cloud (Taurus complex) 
with numerous young stellar objects, while the more turbulent 
(Polaris) is in an inactive complex. 

The exact nature of these intermittent structures, their link 
with shocks, and the role of magnetic fields are still elusive. The 
comparison of observational data with theoretical scalings re- 
quires the ability to compute higher orders of the structure func- 
tions and establish the correspondence between the centroid ve- 
locity and the velocity fields. This stresses the need for large 
homogeneous data samples with at least 10^ spectra. Such data 
sets would also allow determination of the three parameters of 
the class of models to which the SL94 or MHD scalings belong. 

Heterodyne instrumentation {e.g. multi-beam receiver het- 
erodyne arrays) offers a dramatic increase in the spatial dynam- 
ical range accessible, combining high spatial and spectral reso- 
lutions. Sub-arcsecond resolution is needed to resolve the dissi- 
pation scale, combined with a large instantaneous field of view 
to disclose the shape of the dissipative structures. Observational 
signatures of the dissipation of the turbulent kinetic energy might 
be searched for in chemical abundances of species, whose for- 
mation requires high temperatures (Paper I), like CH^, HCO^, 
and water. Excited H2 was also proposed as a good coolant can- 
didate (Falgarone et al. 2005; Appleton et al. 2006). While some 
of these observational requirements are already met by existing 
instruments (e.g. HERA at the lRAM-30m telescope), ALMA, 
SOFIA, and the Herschel satellite will definitely open new per- 
spectives in this field. 
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Appendix A: Two-point statistics 

A.1. Construction of the P„{6C,) 

In each PDF. all the bins which are associated to a number of events less than 
a given value Wmi,,, are blanked. The value of N^\^ depends on the number of 
bins in the histogram. In Fig. lA. 1| we show the P„(dCi) computed for / = 3 in 
the Polaris field, for successive values of A'mi„=0, 10, 30, and 100. It is seen that, 
with A'min=10, the spurious bins having a constant value « 10"* are eliminated. 
The value of each bin and its uncertainty are then determined from the average 
and rms of all the points populating the bin. 



A.2. CVI maps 

The non-averaged CVI map of Fig. lA. 21 computed in the IRAM data for a lag of 
18 pixels (or 180"), shows that large-scale structures exist that are not filtered out 
with large enough lags. The crosses indicate the positions where the CVI values 
of Fig.[T3]have been taken. 

Figure |A3l shows the CVI map computed in the IRAM Polaris and Taurus 
fields, for a lag of 18 pixels. The comparison with the CVI maps of Fig.|9]shows 
that the thin filaments have faded away. However, in the Polaris map, the struc- 
ture visible at a lag of 3 pixels is still visible, though it has broadened. 
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Fig.A.l. VnidCi) computed for Z = 3 in the Polaris field, for 
successive values of Mnm^O, 10 (dashed), 30 (dotted), and 100 
(dot-dashed). 
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Fig. A.2. C VI map computed from the IRAM dataset with a lag 
/ = 18 pixels (or 180"). CVI are in kms"'. The crosses indicate 
the positions used for the cut of Fig. [13] The KOSMA and IRAM 
beams ai^e indicated in the top left corner. 
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Fig. A.3. Same as Fig.|9]for a lag / - 18 pixels. Weak patterns are 
visible, reminiscent of what is seen in Fig|9] Left panel: contours 
are the 0.29 and 0.58 km s"' levels. Right panel: contours are the 
0.29 kms"' with 0.29 kms"' levels. 



IShe et ai] f200lh propose a method of determining tlie parameter f} inde- 
pendently. It makes use of the fundamental assumption of the hierarchical SL94 
model that there is a scaling law for the successive powers of the energy dissipa- 
tion on scale /. Using the Kolmogorov-Oboukhov refined similarity hypothesis 
(RSH, see e.g. iLesieu a recursive relation similar to 

Eq. lA. ll can be written that involves functions of the ratio of successive orders of 
the S „ : 



(A.2) 



with r. 



S p+tiO/S pO). This recursive relation, together with the RSH and 
(oo) ^ leads to the expression of the relative 



the second assumption that J^j 
scaling exponents: 



f;- = pr + (i -3r) 



1 ■ 



1 -q;3 



(A.3) 



with a = and e = 1/3. The a-test proposed bv lSheetal] <2001h is to plot 



r/"+'Vr/-* against 5-,""/^;" since r;"*"/r, 



I 



(Ap/Ai) 



If 



the log-log plot is a line of slope a, the data are said to pass the a-test. In 
Figs. lA.4IA.5l we show the result of the ff-tests applied to the Polaris and Taurus 
data sets. Both pass the a-test with values of a = 0.89 ± 0.01 and 0.950 ± 0.005 
for the Polaris and Taurus fields, respectively. The corresponding values of fS are 
0.70 ± 0.04 and 0.86 ± 0.02 in the Polaris and Taurus fields, respectively (see 
Figs. IA.4IA.5) . The intermittency level in the Polaris field is very close to the 
SL94 value, /3 = 2/3, but our data sets are too small to allow a determination of 
D. 



A.3. Determination of the intermittency level 

I She & Levegud (1994|) developed a model to analyze the small-scale properties 
of an incompressible turbulent flow. Since this model inspired numerous works 
of astrophysical relevance, we summarize its key points here. SL94 propose 
studying the large fluctuations of e/ (defined as the dissipation rate averaged over 



<6f"'>/<6f>. 



balls of size /) through the ratio of its successive moments 

Hence, for each p, the value of ej''' describes the dissipation intensity of a set of 
turbulent structures: as p increases, the associated structures are more coherent 
and more singular. As a result, the hierarchical structures of the SL94 model are 
not related to any physical objects, except for the most intermittent. The SL94 
model has three parameters: the scaling of the velocity with scale /, assumed to 
be that of the K41 theory (v; ~ the level of intermittency characterized 

by a parameter /J; and the dimensionality D of the most intermittent structures. 



Assuming yS = 
to ef"' 



2/3 and D = \, SL94 proposed a recursive relation linking e, 



(P+i) 



(A.l) 



which allowed them to compute the anomalous scaling of the energy dissipation 
rate with scale /, (ep ~ Fp, with Tp = -2/3p + 2[1 - (§/]. 

In principle, by fitting the exponents ((p), it should be possible to determine 
the three parameters of the SL94 class of models. In practice, a reliable determi- 
nation of the three parameters, e.g. by least-square fitting the exponent values, 
requires computing high-order structure functions (p > 6). With 6 orders, we 
could only determine one parameter, the scaling of the velocity fluctuations in 
the cascade v; oc /" and we found 6=1/3 which coincides with the K41 value. 




Fig. A.4. The a-test to both the Polaris field (see Appendix A.3) 
where a - ff with 9 - 1/3. All orders are plotted on the left 
panel and fitted with a single power law, while the right panels 
detail the fits to individual orders p - 2,3,4,5. Fits are done 
over the range indicated by the hne. The individual slopes are 
0.90, 0.91,0.92 and 0.93. 
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Fig.B.l. Cumulative distributions of the CVI {6Ci, full line) and the associated cumulative distributions of 6Cf (dashed line), shown 
as functions of the normahzed 6Ci (computed for a lag / = 3). 
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Fig.A.5. As Fig. IA.4l for the Taurus field, where the slopes for 
each order (/? = 2, 3, 4, 5) are 0.94, 0.95, 0.96 and 0.97. 



Appendix B: Filling factor of the CVI and 
dissipation 

Based on the PDF of CVI of Figs.|4]and|5] we computed the cumulative fraction 
of the CVI SCi and the associated cumulative fraction of SCj. The non-averaged 
CVI SCi are more closely related to the energy dissipation than the azimuthally 
averaged ones ((|(5C;|» because they preserve the two derivatives of the los ve- 
locity. The results are shown in Fig. |B.1| for a lag 1 = 3. Assuming that the 
dissipated energy scales as SCj, we see that in the Polaris field, the regions with 
CVI larger than 3<tsc represent only 2.5% of the surface and they contribute to 
a: 25% of the energy dissipation. In the Taurus field, 5% of the energy dissipation 
is concentrated into less than 1% of the surface. 



